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ABSTRACT 

We present the results of a pilot survey for neutral hydrogen (HI) 21 cm absorption in the Arecibo Legacy Fast 
Arecibo L-Band Feed Array (ALFALFA) Survey. This project is a wide-area "blind" search for HI absorption 
in the local universe, spanning -650 km s~ < cz < 17,500 km s~ and covering 517.0 deg 2 (7% of the full 
ALFALFA survey). The survey is sensitive to HI absorption lines stronger than 7.7 mJy (8983 radio sources) 
and is 90% complete for lines stronger than 11.0 mJy (7296 sources). The total redshift interval sensitive to 
all damped Lya (DLA) systems (N H i > 2 x 10 2l) cm" 2 ) is Az = 7.0 (129 objects, assuming T s = 100 K and 
covering fraction unity); for super-DLAs (Nhi > 2 x 10 21 cm -2 ) it is Az = 128.2 (2353 objects). We re-detect 
the intrinsic HI absorption line in UGC 6081 but detect no intervening absorption line systems. We compute 
a 95% confidence upper limit on the column density frequency distribution function f(NHi,X) spanning four 
orders of magnitude in column density, 10 19 (71/100 K)(l /f) cm -2 < N HI < 10 23 (71/100 K)(l /f) cm -2 , that 
is consistent with previous redshifted optical damped Lya surveys and the aggregate HI 21 cm emission in the 
local universe. The detection rate is in agreement with extant observations. This pilot survey suggests that 
an absorption line search of the complete ALFALFA survey — or any higher redshift, larger bandwidth, or 
more sensitive survey, such as those planned for Square Kilometer Array pathfinders or a low frequency lunar 
array — will either make numerous detections or will set a strong statistical lower limit on the typical spin 
temperature of neutral hydrogen gas. 

Subject headings: radio lines: galaxies — methods: observational — surveys — quasars: absorption lines 



1. INTRODUCTION 

The neutral hydrogen (HI) 21 cm (1420.405752 MHz) spin- 
flip transition has long been employed as an emission line 
tracer of the neutral gas content, surface density profile, kine- 
matics, and dark matter halos of spiral galaxy disks in the 
local universe. HI 21 cm absorption, however, is rarely ob- 
served at low redshift outside of the Galaxy because it relies 
on the chance alignment of a bright radio source with fore- 
ground gas, and the covering fraction of cold neutral gas is 
low. The very slow Einstein rate coefficient for spontaneous 
emission, A2i cm = 2.87 x 10~ 15 s" 1 , means that a detectable 
21 cm line requires either a large gas mass for emission or 
a high column density for absorption; the latter operationally 
limits HI 21 cm absorption line studies to damped systems 
(Nhi > 2 x 10 20 cm" 2 ). But the comoving density of Damped 
Lya (DLA) absorption sys tems is only dN/dz = 0.045 ±0.006 
per unit redshift at z = (IZwaan et al.ll2.Q05h . so identifying 
such systems in the local universe [z < 0.1) requires exami- 
nation of many hundreds of sightlines. If they can be found, 
HI 21 cm absorption lines complement emission lines; while 
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emission line surveys are flux-limited (mass-limited at a fixed 
distance), absorption line surveys are column density-limited 
(at fixed illuminating continuum), independent of distance, 
and can identify very low-mass collections of neutral gas. HI 
21 cm absorption line surveys can potentially identify primor- 
dial atomic gas, infalling or out-going gas in the vicinity of 
galaxies, or spiral disks themselves in a distance- and dust 
extinction-independent fashion. 

The Arecibqj Legacy Fast Arecibo L-Band Feed Array 
(ALFALFA) Survey is a 7000 deg 2 extragalactic HI 21 cm 
emission line sur vey. Compared to the HI Parkes All Sky 
Survey (HIPASS; iMeyer et al .1120041) . the ALFALFA survey 
is 8 times more sensitive, has 4 times the angular resolution, 
3 time s the spectral resolution , and a factor of 1 .4 more band- 
width (Giovanelli et al. 2005), and provides the first data set 
approaching the sensitivity and areal coverage required for a 
truly "blind" HI 21 cm absorption line search. Such a search is 
effectively a search for DLAs that does not rely on UV-bright 
quasar sightlines (and thus can detect dusty DLAs). 

In this paper, we describe the ALFALFA HI 21 cm ab- 
sorption pilot survey, a test of the notion that damped Lya 
systems can be detected "blindly" in the local Universe in 
a shallow, large area spectroscopic survey for HI emission 
from galaxies. We present the methods and results of our 
survey, along with follow-up observations conducted at the 
Green Bank TelescopeQ (GBT), and demonstrate that it is pos- 
sible to place meaningful upper limits on the column density 
frequency distribution function. The pilot survey results sug- 

6 The Arecibo Observatory is part of the National Astronomy and Iono- 
sphere Center, which is operated by Cornell University under a cooperative 
agreement with the National Science Foundation. 

7 The National Radio Astronomy Observatory is a facility of the National 
Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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FIG. 1 . — Radio sources searched for HI absorption in the ALFALFA survey. The blank boxed re gion was not searched due to interference from the strong 
radio source in M 87. The circled source is the re-detected intrinsic HI absorption line in UGC 6081 ( a4.2t . 



gest that a search of the full 7000 deg 2 of the ALFALFA sur- 
vey for HI 21 cm absorption would be worthwhile, both for 
understanding cold neutral gas in the local Universe and in 
preparation for planned absorption line surveys at high and 
low redshift with upcoming new radio telescopic facilities. 

2. SEARCHING FOR ABSORPTION LINES IN THE ALFALFA 
SURVEY 

The ALFALFA survey offers the combined sensitivity, areal 
coverage, and redshift span to make a "blind" HI 21 cm ab- 
sorption line search feasible and astrophysically interesting. 
The complete ALFALFA survey will cover -7000 deg 2 and 
is expecte d to detect more than 2 5,000 galaxies in HI 21 cm 
emission dGiovanelli et al.l 120051) . The data products of the 
survey are dual-polarization 2.4°x2.4° 1024-channel spec- 
tral data cubes spanning 25 MHz each. There are four spec- 
trally overlapping cubes per sky position spanning the range 
-2000 < cz < 17,900 km s" 1 dGiovanelli et al.l 12005). The 
rms noise per spectral chan nel after smoothing to 10 km s" 1 
is 2.2 mJy ( Sainton g^2007l) . 

Our ALFALFA HI absorption pilot survey examines a 517 
deg 2 subset of the mature ALFALFA data product in the 
northern Virgo Cluster region, spanning 10.9 /! < a < 14.95'' 
and +7.7° < 8 < +16.3° (122 spectral cubes). We omit a 1.6 
deg 2 region centered on M 87; the strong radio source in M 87 
(220 Jy at 1 .4 GHz) creates spectral standing waves that ren- 
der this region inaccessible to weak line spectroscopy with 
Arecibo (Giov anelli et al.ll2007l) . 

Figure Q] shows the survey area and distribution of radio 
sources with 8.4 mJy < 5(1.4 GHz) < 5.3 Jy. We searched a 
velocity range -650 km s -1 < cz < 17,500 km s -1 , with no- 
table gaps due to radio frequency interference or Galactic HI 
emission at cz = 0, 2200, 3300, and 8800 km s" 1 of about 
200 km s" 1 width plus a large gap at 15,000-16,000 km s" 1 
dGiovanelli et al.l 120071) . The gaps in the spectral coverage 
render —10% of the velocity range unsearchable. The net 
search covers Az = 0.054 along each line of sight plus the 
range -650 km s -1 < cz < -100 km s -1 . This translates into a 
comoving path length of AX = 0.05710. 

Our line search relies on th e optimized mat ched filter de- 
tection methods developed by Saintonge (2007): we multiply 
the continuum-subtracted data cubes by -1 and run the line- 
finding algorithm as i f we were search ing for narrow Gaussian 
HI emission lines. Saintonge (2007) obtains matched filter 

8 We assume a Qa =0.73 and f2„, = 0.27 cosmology I Komatsu et a l .1201 II) . 
Making no assumptions about global curvature, = l — Q^—Cl m , the co- 
moving path length is AX = Jl= 2 dz(\+z? / \/(\+z) 2 (l+z£lm)-z(z + 2)Cl A . 
We use the AX quantity only in the calculation of the column density fre- 
quency distribution function (Sec. |5J; otherwise, Az is used to describe the 
survey. 



(not single-channel) signal-to-noise (S/N) thresholds for both 
reliability (detections are reliably not spurious) and complete- 
ness (all lines are detected) of 5.5 and 6.5, respectively, for 
narrow lines in the ALFALFA survey. Our absorption line 
search used the conventional ALFALFA matched filter S/N of 
4.6 to search for lines with central depth greater than 7.7 mJy 
toward continuum^ sources with flux densities greater than 
8.4 mJy (this continuum threshold corresponds to a line cen- 
ter optical depth of 2.5 and a neutral hydrogen column density 
of 1.4 x 10 22 cm" 2 assuming hydrogen spin [21 cm line exci- 
tation] temperature T s = 100 K, a source covering fraction / of 
unity, and a 30 km s" 1 FWHM line width). We chose this con- 
tinuum threshold to allow for the possibility of detecting very 
large column densities; while expected to be very rare, the 
number of sources probing such regimes on the sky is large. 
Our low S/N threshold does identify many spurious absorp- 
tion lines, which we subsequently assess using five criteria: 
(1) association with NVSS radio sources (the 45" NVSS beam 
is well-matched to the 1' ALFALFA pixels); (2) similar line 
properties in both linear polarizations; (3) unresolved angu- 
lar size; (4) no association with radio frequency interference; 
and (5) no association with spectral standing waves. Possible 
absorption lines that meet these criteria were subsequently re- 
observed with the GBT (Section©. 

The vast majority of lines detected using our S/N thresh- 
old of 4.6 were obviously spurious, in large part because we 
were searching below the reliability threshold. For the cal- 
culation of column density statistics (Section [5j, we require 
confidence that detected lines are real (reliability) and that a 
non-detection is in fact a non-detection (completeness), so we 
impose the more severe S/N cut of 6.5, corresponding to the 
ISaintongel d2007l) estimate of the 90% completeness limit. For 
the gaussian narrow line template with FWHM of 30 km s -1 , 
this corresponds to a line peak of 11.0 mJy. Allowing for 
opaque lines, we make a continuum cut at the same level. 
For the full line search, there are 8983 radio sources with 
5(1.4 GHz) > 8.4 mJy, but for the restricted column den- 
sity statistics, we reduced the sample to 7296 sources with 
5(1.4 GHz) > 11.0 mJy. The distribution of continua for the 
full sample is shown in Figure|2] 

Supplemental observations were required for two sam- 
ples: candidate HI absorption systems, and strong contin- 
uum sources. Strong sources create spectral standing waves 
between the primary reflector and the superstructure block- 
ing the aperture of Arecibo, frustrating the detection of weak 
lines. It is possible that a technique can someday be devel- 
oped which replicates the "double position-switching" mode 

9 All radio continua used in the selection of sources and the analysis of sur- 
vey results were obtained from the NRAO VLA Sky Sur vey (NVSS) catalog 
and represent integrated flux densities (Condon et al. 1998). 
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FIG. 2. — Flux density distribution at 1.4 GHz of the radio sources searched 
for HI absorption in the ALFALFA Survey. Vertical numbers indicate the 
number of sources in each bin. 

of Arecibo in software using continuum templates within the 
survey data. But in the meantime, it was necessary to observe 
strong continuum sources, 5(1.4 GHz) > 220 mJy, individu- 
ally (Section [3). The strong sources in our sample span the 
range 220 mJy< 5(1.4 GHz) < 5.3 Jy (and exclude M 87 at 
220 Jy), but the majority fall in the range 220-500 mJy. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed 13 candidate absorption systems and 250 
"strong" sources (NVSS flux density greater than 220 mJy) 
with the GBT in April through August, 2008. The GBT is 
particularly well-suited to weak spectral line detection toward 
strong continuum sources thanks to its unblocked aperture. 
Observations of all sources were designed to mimic the spec- 
tral resolution, spectral coverage, and sensitivity of the AL- 
FALFA survey. 

For the follow-up observations of line candidates, we ob- 
served for a total of 10 minutes on-source in a 5-minute 
position-switched mode using two linear polarizations, 9- 
level sampling, 6-second records, and a 50 MHz bandwidth 
centered on the tentative line. The strong sources were ob- 
served in a search mode using two overlapping 50 MHz bands 
spanning 95 MHz from 1335.4 to 1430.4 MHz (cz = -2090 to 
19,080 km s" 1 ) for 5 minutes on-source with a reference "Off" 
spectrum recorded roughly every fifth source. We used two 
linear polarizations, 9-level sampling, and 3-second records. 
For all observations, a calibration diode signal was injected 
for half of each record. 

Each individual off-source-flattened and calibrated spectral 
record was manually examined and flagged for RFI, then av- 
eraged in time and polarization. After polynomial baseline 
subtraction, we achieved 1.3 mJy rms noise in 10.4 km s" 1 
channels (at 1400 MHz) toward the candidate absorbers af- 
ter Hanning and Gaussian smoothing and averaging polar- 
izations. We reached 2.2 ± 0.2 mJy rms noise toward the 
151 strong sources with continua in the range 220-370 mJy, 
2.5 ± 0.3 mJy toward the 83 sources of 370-1000 mJy, 2.8 ± 
0.3 mJy toward the 12 sources of 1 .0-2.0 Jy, and 4.3 ±1.2 mJy 
toward the 4 sources of 2.0-5.3 Jy (the quoted uncertainties 
span the range of rms noise values in each source group- 
ing). All GBT data reduction and analysis was performed 
with GBTIDLB 

4. RESULTS 

10 GBTIDL thttp : //gbtidl ■ nrao ■ eduT) is the data reduction 
package produced by NRAO and written in the IDL language for the reduc- 
tion of GBT data. 
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FIG. 3. — Redshift search path versus HI column density sensitivity. Each 
column density bin indicates the lower bound on the column density de- 
tectable toward the sources in that bin. This distribution represents the re- 
liable and complete sample with a flux density cut of 5(1.4 GHz) > 1 1 mJy 
(as opposed to the 8.4 mJy cut for the sources displayed in Figures \T\ and 
|2). The black, blue, and red histograms indicate the distributions for spin 
temperature modulo covering factor (T s /f) values of 100, 500, and 1000 K, 
respectively. 

4. 1. Column Density Limits 

We made no new detections of HI 21 cm absorption (intrin- 
sic or intervening) in the pilot survey: no ALFALFA absorp- 
tion line candidates were confirmed by GBT observations, 
and no strong sources observed with the GBT show absorp- 
tion lines. Despite the lack of new detections, the uniform 
data set and large number of continuum sources in the sur- 
vey allow us to place limits on the column density frequency 
distribution function f(N,X) (Section|5j. Each non-detection 
of HI absorption toward a radio continuum source places an 
upper limit on the HI column density of 

N(HI) < 1 .8 x 10 18 (T s /f) 1 .064 FWHM • r cm -2 (1) 

where the limit on optical depth is set by S/N= 6.5 correspond- 
ing to the completeness limit, 



r <-ln 1 



11 mJy 
'5(1.4 GHz) 



(2) 



and we assume the minimum FWHM matched filter template 
line width used for the search, 30 km s" 1 . For the strong 
sources observed independently at the GBT, we adjust the 
1 1 mJy line limit above according to the rms noise: it is un- 
changed for the 220-370 mJy sources (2.2 mJy rms), but in- 
creases to 28 mJy for the strongest source at 5.3 Jy. We can 
generally detect any DLA toward the strong sources. Figure 
[5] shows the redshift path searched for each column density 
limit. 

For a source covering fraction of unity and a spin tem- 
perature T, = 100 K, the total redshift interval sensitive to 
all damped Lya (DLA) systems (N H i > 2 x 10 20 cm" 2 ) is 
Az = 7.0 (129 objects); for super-DLAs (N H , > 2 x 10 21 cm" 2 ) 
it is Az = 128.2 (2353 objects). The lowest detectable column 
density under these assumptions is 3 x 10 19 cm" 2 (but see Fig- 
ure 0. 

4.2. Re-detection of Intrinsic HI Absorption in UGC 6081 

We have re-detected the strong intrinsic HI absorption line 
in the interacting system UGC 6081 at cz ~ 10 , 800 km s~\ 
previously detected by iBothun & Schommerl d 19831) and 
IWilliams & Brownl dl983l) . Gaussian fits to the ALFALFA HI 
absorption spectrum (Figure |4|i are listed in Table Q] A two- 
component fit produces a residual spectrum nearly consistent 
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with the rms noise (a three-component fit is not warranted by 
the data). 

UGC 608 1 lies in a region not covered by IRAS, but it does 
appea r in the 2MASS cata log as multiple sources (K s = 12.8- 
13.5. lSkrutskie etaTll2006l) and is likely a luminous infrared 
galaxy (LIRG). Despite the strong HI absorption, the main 
18 cm OH lines were not detected in this system at the level 
of 5.2 mJy rms in 9 km s" 1 channels dKazes & Dickev|[l985l) . 
The FIRST survey resolves this system into two radio sources: 
the NW component is the stronger at 140.3(0.1) mJy; the SE 
component h as an integrated flux density of 30.8(0. 1 ) mJy 
(IWhiteetal.1119971 Figure 13 iBothun & Schommeri (11983k 
and IWilliams & Brownl (1983) assumed that the radio con- 
tinuum was a single source and incorrectly used the aggre- 
gate flux density at 1370 or 1400 MHz to compute the op- 
tical depth and HI column density. Since no interferometric 
spectral line observations have been conducted of this sys- 
tem to date, and the 16" separation of the two continuum 
sources is not resolved by the ALFALFA survey (the ALFA 
beam is 3.3' x 3.8'), the source(s) of the HI absorption remains 
somewhat ambiguous. Note that the measured redshifts of the 
two merging components are consistent with both HI absorp- 
tion l ines and cannot resolve the HI provenance (Geller et al. 
1984). The stronger HI component must originate in the 
stronger (NW) component of UGC 6081 because the line is 
deeper than the flux density of the SE continuum source. The 
weaker, blueshifted HI component, however, could arise in 
either radio continuum source. We compute optical depths 
assuming that all HI absorption originates in the NW radio 
source, and we obtain a very large total HI column density of 
JV(HI) = 1.13(1) x 10 22 (7;/100 K)(l If) cm" 2 by direct inte- 
gration of the observed spectrum (rather than from the Gaus- 
sian fits). This column density is a factor of 1 .4 and 2.8 larger 
than the values obtained by IWilliams & Brownl (119831) and 
IBothun &~S chommeri (119831) . respectively, due in large part 
to larger and confused continua employed in the previous cal- 
culation of optical depths. Note that if T s > 100 K, / < 1, or 
the blue line were associated with the SE component, then the 
already very large column density would be even larger. 

Our re-detection of HI absorption in UGC 608 1 illustrates 
the power of the ALFALFA survey — designed as an HI emis- 
sion line survey — to detect HI absorption (intrinsic and in- 
tervening) in a new completely "blind" survey mode. 

5. ANALYSIS 

We compute the column density frequency d istribution 
function /(Nhi,X) following iCooksevet al. (2010) using the 
limiting column density sensitivity estimates shown in Figure 
® 



f(N m ,X) < 



ANhiAX 



(3) 



where A max is the Poisson upper limit on the detection rate of 
absorption systems with column density Nhi in interval ANhi 
and AX is the comoving absorption path length searched. In 
this case, the search interval is uniform for all sources, so 
AX = 5Xn sens , where SX = 0.0571 is the corrected redshift in- 
terval 5z = 0.0545, and « sens is the number of sources toward 
which observations were sensitive enough to detect a column 
density equal to or greater than Nhi (Figure[3]shows the nearly 
identical quantity Az = &n sen s versus the column density sen- 
sitivity). A 95% confidence upper limit on the Poisson rate is 
Amax = 3.00 when no detections are made. 
Figure [6] shows the column density frequency distribution 
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FIG. 4. — HI absorption spectrum of UGC 6081. The dotted lines show 
the two Gaussian fits to the spectrum, and the upper spectrum shows the 
residual spectrum (data— model) offset by +15 mJy for clarity. The optical 
depth plotted on the right assu mes a continuum flux density of 140 mJy and 
a covering factor of unity (see 34.2t . 





20.0 




10.0 


o 




o 




o 




<N 


10:03:00.0 


3, 




c 




o 






50.0 


a 




c 




o 







40.0 


Q 






30.0 




02:20.0 




Right Ascension (J2000) 

FIG. 5.— Sloan Digital Sky Survey I York et al. 2000) false color image of 
UGC 6081 with FIRST 20 cm radio continuum contours (green; 3.0, 10.4, 
and 36.3 mJy). Red, green, and blue correspond to Sloan z, and r. The 
ALFA beam is larger than the displayed image, but the centroid of the line is 
well-matched to the continuum sources. 

95% confidence limits calculated for a range of T s /f values, 
but in each case a single value is assumed for the entire survey. 
To address this, we also calculated the 95% confidence limits 
on /(Nhi,X) for 1000 trials in which T s /f is uniformly ran- 
domly distributed for each sight line in the range 100-1000 K. 
Figure|7]compares the result of a randomly distributed T s /f to 
the mean of the distribution (550 K), showing little difference 
between the two distributions except in the lowest column 
density bins. This demonstrates that despite the unconstrained 
nature of the spin temperature and covering factor, a distribu- 
tion of values is well-represented in the /(Nhi,X) upper limit 
locus by a mean value. Also, although T s /f is unconstrained 
source-by-source, the prior on this quantity can still produce a 
useful constraint on /(Nhi,X), increasingly so as surveys get 
larger or more sensitive (Section[6]). Figures[6]and[7]compare 
our column density frequency distribution limits to measured 
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TABLE 1 

Intrinsic HI Absorption Lines Toward UGC 608 1 



Line 


V 


V 


z 


FWHM 


Depth 


T 




Nm 




(MHz) 


(kms- 1 ) 




(kms" 1 ) 


(mJy) 




(km s-' ) 


(cm" 2 ) 


1 


1370.943(4) 


10816.3(0.9) 


0.036079(3) 


73(3) 


-46.7(1.6) 


0.40(2) 


31.5(1.6) 


0.57(3) x 10 22 (7/,/100K)(l//) 


2 


1371.14(3) 


10770.8(5.8) 


0.03593(2) 


220(10) 


-17.1(1.4) 


0.13(1) 


30.4(2.7) 


0.55(5) x 10 22 ffi/100 K)(l//) 


Total 


1370.91(3) 


10823.1(5.5) 


0.03610(2) 


86(2) 


-66.5(2.2) 


0.64(3) 


62.1(0.7) 


1.13(1) x 10 22 (7/,/100K)(l//) 



NOTE. — HI absorption properties of UGC 608 1 . We list t he the observed line center frequencies, velocities, and redshifts (barycentric, optical definition), the line widths and 
depths, apparent and integrated optical depths (see Section [4~2l , and HI column density. Values in parentheses indicate uncertainties in the final digit(s). The FWHM and integrated 
optical depth values are in the rest frame of UGC 608 1 . The first two rows in the Table indicate measured and derived quantities obtained from Gaussian fits to the spectrum, while the 
last row in the Table lists numerical values for the entire absorption system obtained without Gaussian fits. 
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FIG. 6. — The column density frequency distribution function, /(Nhi.X). 
All points are 95% confidence upper limits and the horizontal bars indicate 
the column density range spanned by each limit (some arrows and bars are 
omitted for clarity). The long-d ashed line indicates the broken power law 
fit by|Prochaska & Wolfe (2009) fo r absorption systems at z = 2^4, and the 
dotted line indicates a revised fit by Noterdaeme et al. 1 2009). The dash-dot 
line shows the gamma distributi on fit for D LAs in z = galaxies mapped in 
HI 21 cm emission by Zwaan et al. (2005). The black, blue, and red upper 
limits assume spin temperatures modulo covering factor (Ts/f) of 100, 500, 
and 1000 K, respectively. The effect of increasing spin temperature is to shift 
the locus down and to the right because the column density increases linearly 
with T s (Eqn.[TJ. 



f(N HI ,X) distributions from the lZwaan etail (120051) z ~ HI 
21 cm emission line analysis and from the optical survey s of 
redshifted DLAs a t z = 2- 4 by lProchaska & Wolfd (120091) and 
iNoterdaeme et al.l (12009) and show them to be consistent in 
all cases. 

The calculation of f(NHi,X) is influenced by the choice 
of column density interval ANhi because we are calculating 
upper limits: the numerator in Equation [3] is independent of 
the number of sources whereas the denominator grows with 
ANhi- We chose ANhi = 0.5 dex as a compromise between 
good statistics in most bins and resolution of f(Nm,X) along 
the column density axis, but it should be borne in mind that 
these limits do depend on this choice. With coarser sam- 
pling, /(Nhi,X) moves lower but does not become inconsis - 
tent with the fits to f(N H i,X ) ob tained bvlZwaan et all d2005l) . 
iProchaska & Wolfd (120091) . and INoterdaeme et al l (120091) . 

The column density sensitivity is proportional to the combi- 
nation of factors FWHM- T,//; if this quantity decreases, the 
sources are redistributed into lower column bins, driving the 
f(Nni,X) values lower. A line FWHM smaller than the as- 
sumed 30 km s -1 is likely in many DLAs, suggesting that an 
average value of T s /f = 100 K is probably too low (Figure[6]i. 

A source of systematic uncertainty in our calculation of 
f(Nni,X) is the unknown redshift distribution of the illumi- 
nating radio sources. We have assumed for this analysis that 
all radio sources lie at cz > 17,500 km s _I , but there are 
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FIG. 7. — The column density frequency distribution function, f(Nm } X), 
calculated from a uniform random distribution of spin temperatures (mod- 
ulo covering factor) in the range T s /f = 100-1000 K. All points are 95% 
confidence upper limits and the horizontal bars indicate the column density 
range spanned by each limit (some arrows and bars are omitted for clarity). 
The vertical error bars indicate the 1 a range in the 95% confidence upper 
limit obtained from 1000 trials. For reference, the blue line indicates the 
upper limit obtained for a uniform spin temperature of 550 K (the mean of 
the distributio n). The long-dashed line indicates the broken power law fit by 
IProchask a & Wolfe ( 2009) for absorption systems at z = 2^1, and the dotted 
line indicates a revised fit by [Noterdaeme et al. 12009). The dash-dot line 
shows the gamma distribu tion fit for D LAs in z = galaxies mapped in HI 
21 cm emission by Zwaan et al. ( 2005). 

certainly cases to the contrary (e.g., UGC 6081). This er- 
ror operates in one direction: it diminishes the total redshift 
search path of the survey, driving our f(Nm,X) limit upward. 
The size of this effect, however, is very small: using the red- 
shift distributio n of NVSS sources w ith 5(1 .4 G Hz) > lOmJy 
determ ined by iBrookes et al.l d2008l) and fit by Ide Zotti et alj 
d2Ol0h . we calculate that no more than about 60 sources in the 
517 deg 2 pilot survey lie within the observed redshift range, 
which is 0.8% of the sample. The redshift distribution of these 
sources will also be skewed toward the high redshift side of 
the survey, so the effect of this incomplete redshift coverage 
is truly negligible. 

Since we obtain only upper limits on the column density 
frequency distribution /(Nhi,X), its moments, such as the HI 
mass density (pui or equivalently ilm) and the comoving cov- 
ering fraction £(X), are dominated by the endpoint limits and 
thus not well constrained by our observations. It is likely that 
such moments can be computed from future surveys (Section 
16.11 ). Our lack of detections is perfectly consistent with expec- 
tations and recent work on damped Lya lines and HI 21 cm 
emission (Section|6]). 

6. DISCUSSION 

Our upper limits on the column density frequency dis- 
tribution function are consistent with previous determina- 
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tions of f(N H i,X), both at hig h redshift dProchaska & Wolfd 
2009; Noterdaeme et al. 2009, via Lya absorption) and low 
dZwaan et al.ll2.Q05l via HI 21 cm emission). Figures [6] and [7] 
show that an expanded survey would either start making ab- 
sorption line detections or will place a strong statistical limit 
on the hydrogen spin temperature to source covering fraction 
ratio (T s /f). We suspect that it will be the latter, but in any 
case, our results provide strong motivation for a more sensi- 
tive or a larger area survey. For example, a search for ab- 
sorption lines in the full ALFALFA survey would increase the 
comoving search path AX by a factor of 14 (1.15 dex) in each 
Nhi bin, which will place the f(Nm,X) 95% confidence lim- 
its in conflict with previous measurements of the distribution 
over the range 10 20 cm" 2 < Nhi ^ 10 22 cm" 2 unless either 
DLAs are detected or T s /f > 500 K on average. 

Clearly a more productive approach to obtain the column 
density frequency distribution f (Nhi, X) at z = is via HI 
21 cm emission lines (e.g., Zwa an et al.H2.Q05h . but beyond 
0.2, where the emission line flux density becomes pro- 
hibitively weak for anything less than a square kilometer of 
collecting area, one would use a distance-independent absorp- 
tion line survey such as the pilot survey presented here. Fur- 
thermore, "blind" 21 cm-based detections of DLAs can iden- 
tify dusty sightlines that would otherwise be missed by opti- 
cal/UV surveys, which is key to identifying exceptionally rare 
molecular absorption systems and addressing the fraction of 
dusty DLAs. 

6.1. Future Surveys 

As Figures [3] and [7] demonstrate, even incremental im- 
provements in survey sensitivity would lead to large improve- 
ments in column density sensitivity and statistics because the 
source population grows as a power law as the 1 .4 GHz flux 
density decreases (Figure |2). But for strong sources, which 
provide probes of low column density, areal coverage remains 
important. 

The first step should be to search the full ALFALFA survey 
for 21 cm absorption lines; this pilot survey covered 7% of 
the full survey, so an order of magnitude increase in search 
path is possible using extant data. This will either produce 
detections of intervening absorption line systems or lower the 
limits on the column density frequency distribution function, 
indicating a statistical lower bound on spin temperature (mod- 
ulo covering fraction). It is already clear from Figure [6] that 
the typical spin temperature of HI gas at z ~ cannot be less 
than T s /f ~ 100 K. 

Future facilities that could perform commensal or dedicated 
21 cm absorption line surveys include Square Kilometer Ar- 
ray (SKA) prototypes, including the 256-element Allen Tele- 
scope Array, the SKA Molonglo Prototype, the Australian 
SKA Pathfinder, and the Karoo Array Telescope (MeerKAT), 
as well as Epoch of Reionization (EoR) and Dark Ages tele- 
scopes, such as a lunar far-side array, the Low Frequency 
Array (LOFAR), the Precision Array to Probe the EoR, and 
the Murchison Widefield Array. These facilities will provide 
various combinations of areal coverage, bandwidth (Sz per 
source), and sensitivity. As we have demonstrated here, sen- 
sitivity is key because it determines the number of sightlines 
available per solid angle of sky surveyed, and the power-law 
increase in radio source counts with decreasing flux density 
means that even slight improvements in sensitivity will lead to 
large gains in the size of absorption line surveys. Bandwidth 
and areal coverage, in comparison, grow absorption line sur- 



veys linearly. 

While absorption line detection is distance independent, 
which is an enormous advantage over redshifted emission line 
surveys, future HI absorption line surveys will have to con- 
tend with the issue of unknown redshifts of the continuum 
population and the loss of foreground sources in the total red- 
shift search path. 

7. CONCLUSIONS 

Our "blind" pilot survey for HI 21 cm absorption lines 
in a 517 deg 2 section of the ALFALFA survey spanning 
-650 km s" 1 < cz < 17,500 km s" 1 detected no new absorp- 
tion line systems, but re-detected the intrinsic absorption in 
UGC 6081. We calculate column density sensitivity limits 
toward a sample of 7296 radio sources and construct up- 
per limit envelopes for the column density frequency dis- 
tribution function f(Nm,X) that are co nsistent with previ- 
ous work at high and low redshift (e.g.. IZwaan et al] 120051 : 
iProchaska & Wolfd l2009t iNoterdaeme et alJl2009l) .~Our re- 
sults suggest that higher redshift, larger solid angle, larger 
bandwidth, and especially more sensitive surveys will be an 
astrophysically worthwhile pursuit in terms of constraining 
the 21 cm spin temperature and detecting dusty damped Lya 
systems as well as DLAs toward optically faint but radio- 
loud illumination sources. A subset of these HI 21 cm ab- 
sorption systems will also host the exceptionally rare and elu- 
sive molecular absorption systems, which can be employed to 
study cold molecular gas at high redshift and to constrain the 
cosmic evolution of fundamental physical constants. 
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